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ABSTRACT 

This paper presents numerical simulations of test particle Fermi acceleration at relativistic 
shocks of Lorentz factor r sn = 2 — 60, using a realistic downstream magnetic structure 
obtained from the shock jump conditions. The upstream magnetic field is described as pure 
Kolmogorov turbulence; the corresponding downstream magnetic field lies predominantly in 
the plane tangential to the shock surface and the coherence length is smaller along the shock 
normal than in the tangential plane. Acceleration is nonetheless efficient and leads to power- 
law spectra with index ~ 2.6 — 2.7 at large shock Lorentz factor T s h 1, markedly steeper 
than for isotropic scattering downstream. The acceleration timescale i acc in the upstream rest 
frame becomes a fraction of Larmor time in the ultra-relativistic limit, i acc rs lOt^/T^. 
Astrophysical applications are discussed, in particular the acceleration in 7— ray bursts inter- 
nal and external shocks. 
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1 INTRODUCTION 

The Fermi acceleration process of charged particles bouncing back 
and forth across a shock wave is the main ingredient for the gener- 
ation of high energy radiation in a variety of astrophysical environ- 
ments. This observed radiation is generally synchrotron light emit- 
ted by the accelerated electrons; in this case one may recover the 
spectral index s of the accelerated population from the synchrotron 
index. For example, the afterglow emission of 7— ray bursts that is 
seen in X-ray through the infrared is generally interpreted as syn- 
chrotron emission of electrons accelerated at the ultra-relativistic 
external shock of Lorentz factor T s h ~ 300. The inferred spectral 
index, s ~ 2.3 ± 0.1 (Waxman 1997; see also Meszaros 2002 and 
Piran 2004 for reviews), thus probes the nature of shock accelera- 
tion in the ultra-relativistic regime. Similarly synchrotron emission 
of electrons accelerated in the mildly relativistic internal shocks 
(r s h ~ 2 — 5 in the comoving frame) with index s ~ 2.3 ± 0.1 
could explain the prompt 7 emission (see e.g. Meszaros 2002, Pi- 
ran 2004 and references therein). These observations thus provide 
anchor points for studies of Fermi acceleration in the moderate to 
the ultra- relativistic regime. They have actually been regarded as 
a dramatic confirmation of the theory of shock acceleration in the 
relativistic regime, which has been claimed to predict a "universal" 
asymptotic spectral index s ~ 2.23 in the ultra-relativistic regime 
r sh > 1. 

Relativistic shock acceleration has been studied through a va- 
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riety of methods, either analytical (Peacock 1981, and more re- 
cently Vietri 2002, Vietri 2003, Keshet & Waxman 2005, Blasi & 
Vietri 2005), semi-analytical (Kirk & Schneider 1987; Gallant & 
Achterberg 1999; Kirk et al. 2000; Achterberg et al. 2001), or nu- 
merical (Ellison et al. 1990, Ostrowski 1991, Ballard & Heavens 
1992, Ostrowski 1993, Bednarz & Ostrowski 1996, 1998, 1999, 
Ellison & Double 2002, 2004, Lemoine & Pelletier 2003, Meli & 
Quenby 2003a,b, Bednarz 2004, Niemec & Ostrowski 2004, Bar- 
ing 2004). Not all of these studies find the universal value for s, 
however, all the more so when anisotropic configurations such as 
oblique shocks are considered. One clear example is the demonstra- 
tion that Fermi acceleration in superluminal (perpendicular) shocks 
in the absence of cross-field diffusion becomes inefficient (Begel- 
man & Kirk 1990); in the relativistic regime, oblique shocks are 
superluminal unless the angle between the magnetic field and the 
shock normal Ob < l/r s h. 

It is generally suspected that the inclusion of scattering would 
make Fermi acceleration more efficient in the relativistic regime. 
The simulations of Bednarz & Ostrowski (1998) and Baring (2004) 
have indeed confirmed that the spectral slope tends to increase with 
increasing shock obliquity and with decreasing turbulence level, 
whose role is to permit cross-field line transport to the shock front. 
However, a limitation of these simulations is that the scattering is 
simulated in a phenomenological way by setting a ratio of the per- 
pendicular to parallel diffusion lengths and drawing pitch angles 
at random at each time step. Upstream, it has been demonstrated 
that an ultra-relativistic shock wave overtakes the particle before 
this latter has had time to scatter efficiently (Gallant & Achterberg 
1999, Achterberg et al. 2001), so that the details of particle trans- 
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port are probably not crucial. Downstream, however, the particle 
has to turn back before re-crossing the shock, and the approxima- 
tion of ad-hoc diffusion lengths may be too naive to accurately sim- 
ulate the transport. In fact, one may expect non-trivial correlation 
functions between displacements along different directions as well 
as subdiffusion regimes to play a significant role in the return to the 
shock. 

Several studies have tried to integrate out exactly the particles 
trajectories in a well-defined magnetic field structure (Ballard & 
Heavens 1992, Ostrowski 1993, Lemoine & Pelletier 2003, Niemec 
& Ostrowski 2004). The study of Ballard & Heavens (1992) in- 
volved a realistic magnetic field structure, in the sense that it obeys 
the shock jump conditions, but was limited to mildly relativistic 
shocks (r s h ^ 5). Nevertheless, it observed a trend of increas- 
ing spectral index with increasing shock velocity, a result which 
has been disputed by the more exhaustive simulations of Ostrowski 
(1993). The recent work of Niemec & Ostrowski (2004) consid- 
ered situations of moderate turbulence levels with varying degrees 
of obliquity for mildly relativistic shocks {Tsi\ ^ 5). The con- 
clusions obtained indicate that various spectral slopes can indeed 
be obtained, although the noise on the simulations is not negligi- 
ble. Finally, the work of Lemoine & Pelletier (2003) introduced 
a new numerical Monte Carlo method to study relativistic Fermi 
acceleration, on which the present work is based. It assumed the 
downstream turbulence to be isotropic and confirmed the value 
s ~ 2.2 — 2.3 predicted in that case. 

It seems fair to say that a clear picture of the efficiency of 
relativistic Fermi acceleration in a magnetic structure that includes 
compression of the upstream magnetized configuration has not yet 
emerged. The present paper proposes to undertake such simulations 
in order to make progress along these lines. We assume that the up- 
stream magnetic field is described by pure Kolmogorov turbulence, 
i.e. there is not uniform component; this can be seen as the limit 
SB / B — > oo of a highly turbulent plasma. The upstream magnetic 
field is compressed by the shock into an anisotropic downstream 
turbulence. We conduct our simulations in the mildly and ultra- 
relativistic regimes; simulations of that kind in this latter regime 
have never been attempted before. 

In Section 2, we describe in detail the numerical techniques 
and the procedure used (borrowed from Lemoine & Pelletier 2003) 
to simulate the Fermi acceleration process. In Section 3, we present 
our results on the (downstream) return probability, the accelera- 
tion timescale and the accelerated spectrum as a function of shock 
Lorentz factor. In Section 4 we discuss the relaxation length of the 
turbulence and argue that, for relativistic shocks at least, particles 
that return to the shock downstream do not travel beyond the point 
where the anisotropy of turbulence has relaxed. We also discuss 
the properties of transport of particles in the strongly anisotropic 
turbulence generated by shock compression, compare our results 
to previous studies and comment on the applications of our results 
to shock acceleration in 7— ray bursts and to shock acceleration of 
ultra-high energy cosmic rays. Conclusions and a summary of the 
results are provided in Section 5. 



2 NUMERICAL SIMULATIONS 

2.1 Jump conditions and magnetic fields 

In the present work, we assume that the magnetic field is dynami- 
cally unimportant, i.e. its energy density can be neglected with re- 
spect to that of the fluid. We also consider a strong shock, for which 



the upstream random kinetic energy per particle can be neglected 
with respect to that downstream. The corresponding hydrodynamic 
jump conditions are given in Blandford & McKee (1977), and re- 
viewed in Kirk & Duffy (2001) and Gallant (2002). The shock 
Lorentz factor is denoted T s h in the upstream frame (taken as the 
lab frame), and the shock velocity upstream is /3 s h. Unless other- 
wise noted, all quantities are calculated in this frame. If relevant the 
reference frame is indicated by a subscript, e.g., /3 a h| d refers to the 
shock velocity measured in the downstream rest frame and r sh | d 
refers to the shock Lorentz factor in the downstream frame. The 
relative Lorentz factor between upstream and downstream is noted 
r ro i and reads: 

r ro i = r S hr sh | d (i — /3 S h/3 S h|d)- 

The downstream Lorentz factor r sh | d as well as r rc i can be ob- 
tained as a function of V s h (upstream shock Lorentz factor) using 
the relations derived from the shock jump conditions for a Synge 
equation of state (Gallant 2002): 
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where £ = mc 2 /Ta, T d being the downstream temperature and 
m the particle mass, G(£) = Kz (£) /A"2(£)> wrtn K2, K3 modi- 
fied Bessel functions, and = [£(?(£) - l] 2 - £ 2 . These rela- 
tions hold for a gas composed of possibly different particles species 
but with same £ (Gallant 2002). Equations Q 

can be inverted nu- 
merically to obtain F sll | d as a function of r s h. In particular, in the 
ultra-relativistic limit T s h — * +00, one finds the well-known re- 
sults /3 sh|d -> 1/3 (r sh]d -> 3/V8) and T rc i F sh /V2. 

The conservation of the electromagnetic field energy- 
momentum tensor implies the following jump conditions for the 
magnetic field components By (aligned with the shock normal) and 
B± (tangential to the shock surface): 
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and as before, /3 sh | u = /3 s h, r sh | u = T s h. The parallel com- 
ponent B11 is thus conserved while the perpendicular component 
B± is amplified by the proper shock compression ratio R = 
/3sh| u r s h|u/,3sh|dr s h| d . In the ultra-relativistic limit R — > r sh v^, 
and the total magnetic field strength is amplified by W2/3R. 

We assume that the upstream magnetic field is purely turbu- 
lent with a power spectrum describing Kolmogorov turbulence with 
maximal length scale L max . It is modeled as a sum of static plane 
wave modes according to: 



Bu(x) = 



(3) 



with ek a unit polarization vector orthogonal to k, <jik a ran- 
dom phase and |Gfc| 2 oc k~ 5 ^ 3 the amplitude of the power 
spectrum. The wavenumbers k range from fc m i n = 27r/L max to 
some maximal wavenumber fc max ^> fc m i n ; numerically we em- 
ploy 250 wavenumbers modes whose directions are drawn at ran- 
dom, and whose moduli are spaced logarithmically between fc m i n 
and fc max = 5 • lO^fcmm- The amplitude Gk can be chosen as real 
and is normalized such that: 



I / dxB 2 (x) = £|G fc | 2 = B 



2 
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(4) 



with B lms the squared turbulent magnetic field strength. 



According to the shock jump conditions, the downstream 
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magnetic field is described by an anisotropic turbulence: while Sii 
is conserved, the turbulence wavenumbers fcii are amplified by R, 
which corresponds to the compression of the eddies by 1/R along 
the shock normal. The perpendicular wavenumbers k± are con- 
served through the shock but B± is amplified as before. Hence the 
downstream magnetic field is described by: 



Bd(x) = £ 



(5) 



where k is related to the wavenumber k of Eq.|5|by ky = ^ky 
and kx = kx; similarly §k is related to ek by: eiit = en k and 
ex,k = f?ex,k- Note that k ■ e k = k ■ ek = as required for a 
divergenceless field; 4>^ and Gk are not modified. We chose not to 
normalize the above polarization vector to unity downstream and its 
modulus gives the overall amplification factor of the magnetic field. 
This is but a matter of convention: one may equally well embody 
the compression factor in Gk and normalize e k to unity. 



2.2 Monte Carlo simulations 

Our numerical procedure is summarized in Lemoine & Pelletier 
(2003). It consists in two main steps: in a first stage, we conduct 
Monte Carlo simulations of particle propagation in a magnetized 
medium (either upstream or downstream) and derive the statistical 
properties related to shock crossing and re-crossing, as described 
below. In a second step we use these statistical distributions in con- 
junction with the Lorentz transforms from one frame to the other to 
reconstruct the accelerated spectrum that escapes downstream. 

Once the magnetic field structure is set up as described above, 
one Monte Carlo simulation of the propagation of particle con- 
sists in integrating the equation of motion in the magnetic field. 
The particle trajectory is saved in time intervals that are a fraction 
f u ~ 10 -4 (upstream) or /d — 10~ 2 (downstream) of Larmor time 
£l = Rl/c (with RrmL = p/qB) over a time period as long as 
AT U ~ 10 2 (upstream) or ATd ~ 10 4 (downstream) Larmor time. 
For each computed trajectory one can build a statistical sample of 
shock crossing and re-crossing as follows. One draws at random a 
point along the trajectory, which defines the point at which the par- 
ticle enters through the shock. One records the ingress pitch angle 
cosine of the particle momentum with respect to the shock normal 
at that point. One then searches for the point along the trajectory 
at which the particle exits through the shock and the corresponding 
egress pitch angle cosine is recorded. In the downstream medium, 
it happens that the particle never re-crosses the shock as the shock 
itself moves away with speed /3 sh | d ~ 1/3 (r s h S> 1). By iterating 
the above procedure, i.e. drawing other points of entry in the trajec- 
tory, and building other trajectories, one can measure the probabil- 
ity laws that control in a direct manner the Fermi process. 

In particular, the ratio of the number of shock re-crossings to 
the total number of shock entries at a given ingress "pitch" an- 
gle (defined here as the angle between the momentum and the di- 
rection of the shock normal) cosine fi' gives the return probabil- 
ity P le t(p')- In a similar way, the number of shock re-crossings 
through an egress pitch angle cosine /i c for a given ingress cosine 
p 1 gives (after proper normalization) the conditional return proba- 
bility V(p'; p c ). One can define and calculate these quantities both 
downstream, Vd(pd\ Pd)' an d upstream, V u (p u ; p u ). Note that the 
ingress and egress pitch angles are calculated in the rest frame of 
the fluid under consideration. The normalization of the conditional 
probability laws is such that their sum over the egress pitch angle 
cosine yields the return probability as a function of ingress pitch 



angle cosine: 



fret,dO"d) = J dp C d'Pd(Pd;Pd), 

Pict,u(Pu) = / dPuVuiPu^Pu)- 



(6) 



Obviously the upstream return probability P rc t,u must be 
unity if one considers an infinite planar shock with an infinite life- 
time. This provides a useful check on the numerical procedure; in 
the present calculations, P Ic t. u does not deviate from unity by more 
than ~ 10 -6 . Downstream it is mandatory to verify that one does 
not miss possible late returns by varying the trajectory integration 
time; we estimate that the mean of P rc t,d over ingress pitch angle 
cosines is accurate to better than ~ 10~ 4 . 

Finally these simulations give a direct measurement of the re- 
turn timescale to the shock as a function of pitch angles. This mea- 
surement is particularly important to estimate the maximal accel- 
eration energy in a variety of environments, as discussed in Sec- 
tion 3.3. 

Once the upstream and downstream laws of return probability 
are known, the simulation of the acceleration process itself can be 
performed as follows. We denote by J r ^ n+1 (/Xd, £d) the distribu- 
tion function of particles that enter the shock towards downstream 
with ingress pitch angle cosine /id, that have experienced 2n + 1 
shock crossings and that carry energy e<a (downstream frame). Sim- 
ilarly we define the distribution function jT 2n (/i u , e u ) of upstream- 
going particles with ingress pitch angle cosine /i u , having experi- 
enced 2n shock crossings and carrying energy e u . If we denote by 
T° the injection population upstream, then after an even (resp. odd) 
number of shock crossings the particles are necessarily upstream 
(resp. downstream). The injection (isotropic) distribution function 
jFu is normalized to unity, as follows: 



d/i u de u 



= 1. 



(7) 



The integral over /i and e of the distribution functions J r ^ n+1 (ji, e) 
and jT 2n (/i, e) with n > is smaller than unity, due to escape of 
particles downstream at each cycle. 

Now, particles that enter upstream after 2n shock crossings 
with ingress cosine /i u re-cross the shock with egress cosine p^ 
and with conditional probability V u (p u > Pu)- The total number of 
particles with egress pitch angle p 1 ^ and energy e u at the 2n + 1 
shock crossing is f dfj^ V u (p u ', Pu)3~u n (pu, e u ). We note that the 
upstream egress cosine p^ is related to the corresponding down- 
stream ingress cosine [i d , by a Lorentz transform, just as the en- 
ergies measured in the upstream frame (e u ) or downstream frame 
(ea): 



Pd 



pu — Prcl 
1 - A-cIMS : 



£d — r re i(l — /3 le iAt u )e 



(8) 



with a similar relation between p d and ^ when the particle crosses 
the shock from downstream to upstream. 

Therefore, the conservation of particle number at shock cross- 
ing u^d implies the following relation between J r ^ n+1 and T^ 1 : 



d" +1 (Md, £d) d/i d de d 



dp u V u (p u ; pD^u" {Pu, e u ) 
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d^°de u , 



(9) 



and one obtains a similar system for shock crossing d— >u: 
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with: 



Md + /3rcl 
1 + /3rel/i d ' 



= r re i(l+/3 relA t d )e d (11) 



where the " ~ " symbol has been introduced to differentiate the 
values of /i d and e d from one cycle (2n — 1 shock crossings) to 
the next (2n + 1 shock crossings). The integration bounds on p are 
imposed by the shock crossing conditions. 

The terms within brackets in Eqs. {5) and i lOt correspond to 
the distributions upon exit from upstream and downstream respec- 
tively. These equations assume implicitly that the conditional prob- 
ability laws do not depend on energy. This will be shown to be a 
good approximation in Section 3.1. 

After each cycle u^d^u, a population jF " t +1 (e d ) = 
/ d/i d [l — P ro t(A'd)]-7 : d Tl+1 (Md; £ d) of the particle population has 
escaped downstream. The sum over n of these escaping particles 
forms the outgoing accelerated particle population: 



^out(e) = ^ ■7 r out +1 (e) 



(12) 



By following each shock crossing, and using Eqs. ©, OD, (TO), 
m and o one can follow the evolution of J~d, J~u and J~out, 
starting from a mono-energetic and isotropic initial injection distri- 
bution upstream. A similar formal development of the acceleration 
process by repeated shock crossings has also been proposed inde- 
pendently by Vietri (2002): the flux of particles crossing the shock 
in the stationary regime, noted Ji n in Vietri (2002) is related to the 
above as J ln = CJ2™~Zq°° T^ a+1 with C a normalization con- 
stant (see also Lemoine & Pelletier 2003). 

The present technique has significant advantages when com- 
pared to standard Monte Carlo techniques which follow the particle 
trajectories on both sides of the shock through the whole acceler- 
ation process; in particular, it offers a significant gain in signal to 
noise as will be obvious in Section 3. 

It has however one caveat that should be underlined and which 
resides in the fact that we compute the accelerated spectrum by 
merging separate pieces of information on transport properties up- 
stream and downstream. In so doing, we neglect the possible cor- 
relations that may exist between the upstream magnetic configura- 
tion at the crossing point and that downstream, i.e. we neglect the 
possible correlations between upstream and downstream transport. 
The only method that can take this effect into account is the di- 
rect Monte Carlo integration of individual particle trajectories. It 
is therefore important to compare the results obtained with these 
two methods in order to assess the magnitude of this effect. In 
the case where scattering is isotropic downstream, it appears that 
various methods converge to the same value of the spectral index, 
and this includes various non-Monte Carlo methods which cannot 
take the above effect into account (see, e.g., Achterberg et at 2001, 
Lemoine & Pelletier 2003, Keshet & Waxman 2005) as well as di- 
rect Monte Carlo methods (e.g., Bednarz & Ostrowski 1998). This 
suggests that, at least in the isotropic limit, these correlations do not 
play a significant role in the determination of the spectral index. 
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Figure 1. Downstream return probability vs ingress pitch angle cosine /Li*, 
(downstream rest frame) for various shock Lorentz factors, as indicated, and 
for a rigidity 2n i? L /L m ax = 6 • 10~ 4 . 



3 RESULTS 

The numerical technique described in the previous section allows 
to collect a significant amount of information on the acceleration 
process, in particular the conditional probabilities of return from 
downstream or upstream, the energy gain per cycle as well as the 
acceleration timescale. In order to better understand the results ob- 
tained for each of these quantities, it is necessary to emphasize the 
difference between the effective coherence length along the shock 
normal L« and that tangential to the shock front L± as measured 
downstream, see Section 2: Lu — L±/R = L mBX /R, where L max 
is the coherence length of the upstream magnetic field, and R the 
proper shock compression ratio, R ~ r s h\/8 when L s h 3> 1. 
This distinction takes on a particular importance when one com- 
pares the results over various values of the shock Lorentz factor 
and over various values of the rigidity p = 2ivRi J /L lna , x , where 
i?L denotes the Larmor radius in the rest frame of consideration. 
In principle, the transport properties of particles in a magnetic field 
depend solely on the rigidity. However, when the effective coher- 
ence length along the shock normal depends on L s h through R, 
while L± does not, there is no unambiguous definition of rigidity. 
In particular, the above definition of p does not correspond to the 
effective rigidity pu = 2nRi,/Lu that controls the scattering of 
particles with turbulence modes of wavevector parallel to the shock 
normal: a given rigidity p = 2-kRl/ X max corresponds in fact to 
larger and larger values of p« as T s h increases. The relevance of 
this observation to the results will be adressed shortly. 



3.1 Return probability 

The average return probability marginalized over egress angle, 
Pret,d> defined in Eq. [6| as the direct average of the conditional 
probability law Vd(pd, Pd) over me egress pitch angle cosine /i d , 
is shown as a function of the ingress pitch angle cosine p d in Fig.Q 
The increase of P TCt ,d as p d — > 1/3 backs up the notion that parti- 
cles crossing the shock from upstream to downstream at near graz- 
ing incidence with the shock front have a substantially higher prob- 
ability of returning to the shock than those crossing the shock head- 
on (p d -> -I)- 

When viewed as a function of shock Lorentz factor, the av- 
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Figure 2. Downstream return probability averaged over ingress angle vs 
shock Lorentz factors for various rigidities as indicated. The dotted line is 
an ad-hoc fit (P ro t) = 0.97 — 0.67/3 s h which accounts well for the depen- 
dence of the return probability with shock Lorentz factor at low rigidities. 



erage return probability appears to reach an asymptotic law, as is 
made apparent in Fig.Q This is not a trivial result in itself, as the 
nature of the turbulence downstream depends rather strongly on the 
shock Lorentz factor. Section 4.1 provides examples of downstream 
trajectories for two different values of F s h, and indeed, the dis- 
placements along the shock normal differ widely. Hence one might 
naturally expect that the return probability would carry some form 
of dependence on the shock Lorentz factor r s h. As we now argue, 
this is related to the fact that the scattering timescale in the direction 
along the shock normal, i.e. the time required for the particle to turn 
back, is a function of Larmor time, as demonstrated in Section 4.1. 
There it is argued that the particles that return to the shock have 
done one reflection on the compressed turbulence in their first in- 
teraction; indeed, particles get trapped in a layer of the compressed 
turbulence when r s h ^> 1 and p <C 1, hence they cannot return to 
the shock unless they do so in the first interaction. This reflection 
is in fact a half-gyration of the particle around a field line which 
is mainly oriented along the shock front as a result of shock com- 
pression, and this explains why the scattering time is of order of 
the Larmor time. For large values of r s h, the shock velocity with 
respect to downstream (3 Bh ^ d -» 1/3 becomes independent of r s h, 
and so does the scattering timescale (at a given rigidity, see Section 
4.1), hence so does the return probability. 

One may further average the return probability P re t,d over the 
ingress pitch angle in order to define the average return probability 
fret : 



(Pret) 



/d^P rct ,d(Md) 



(13) 



This probability is shown as a function of shock Lorentz factor for 
varying values of the rigidity p in Fig.|2| The dotted line represents 
the empirical fit: (P ro t) — 0.97 — 0.66/3 srl which provides a good 
approximation at low rigidities. This figure shows how the average 
return probability reaches an asymptote with r s h for sufficiently 
low rigidities. At high rigidities (upper panels), one recovers a de- 
pendence of (P re t) on F s h. This latter effect is likely related to the 
factor R difference between p« and p: as pii = Rp becomes larger 
than ~ 0.1 — 1, particles can no longer interact resonantly with the 
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Figure 3. Downstream return probability averaged over ingress angle vj 
rigidity for various shock Lorentz factors, as indicated. 



turbulence wave modes (Casse et al. 2002); they take a longer time 
to return to the shock, see Section 3.3, and their return probability 
becomes sensitive to the nature of the turbulence, hence to r a h. 

Finally one can plot the average return probability as a func- 
tion of rigidity for various values of the shock Lorentz factor, see 
Fig. It is important to note that the average return probability 
does not depend on the rigidity, at least for sufficiently low rigidi- 
ties p« <g 1 for the same reasons as above. The conditional return 
probabilities (from which (P rc t} is obtained) are also found not 
to depend on rigidity in that range. In order to measure the spec- 
tral index of the accelerated spectrum for rigidities in the inertial 
range of resonance, it is important to use only the datasets of the 
smallest rigidities downstream, i.e. p = 6 ■ 10 -4 and p = 2 • 10 -3 , 
where we still have fc max PL > 1, for the same reasons as discussed 
above. For upstream probability laws, one can use all datasets since 
the rigidities are well in the inertial range in the absence of com- 
pression effects. In what follows, we use different combinations 
of one downstream with one upstream of these datasets to simu- 
late the Fermi acceleration process and measure the spectral index. 
We use these different datasets as independent realizations of the 
conditional probability laws in order to estimate the numerical un- 
certainty on the spectral index. 



3.2 Escaping accelerated particles 

The fraction of particles that do not return to the shock adds up to 
form the outgoing accelerated particle spectrum. As shown by Bell 
(1978) for the case of non-relativistic shocks, the spectral index 
of this spectrum is determined by the average return probability 
and the mean energy gain at each cycle. For relativistic shocks, 
the analytical development of Vietri (2002), whose formulation is 
very similar to that presented in Section 2, shows that the spectral 
index is determined by the energy gain properly averaged over the 
equilibrium distribution functions both upstream and downstream, 
see also Lemoine & Pelletier (2003). 

The average energy gains per cycle u^d^u and half-cycles 
d^u, u-^d are shown in Fig. [4] which shows clearly that the gain 
is of order ~ T^ h for the first complete cycle u— >d— >u, and falls 
to < 2 in subsequent cycles, as anticipated by Gallant & Achter- 
berg (1999), Achterberg et al. (2001). This strong limitation of 
the energy gain is due to the anisotropy of the distribution func- 
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Figure 4. Average energy gain per cycle u— >d— >u (diamonds), per half cy- 
cle u— >d (triangles), and per half cycle d^u (squares) plotted vs successive 
cycles. 



tion upstream: particles do not have time to be deflected by an an- 
gle greater than ~ 1 /T s h upstream before being overtaken by the 
shock which moves at speed /3 s h — 1 with respect to upstream, 
hence the particles energy is decreased by a factor l/r re i in the 
half-cycle u^d through the Lorentz transform, in agreement with 
Fig.|l| In order to return to the shock downstream, particles must 
turn back and the average energy gain in the half-cycle d^u is 
now ~ «r re i, with k < 2 resulting in the total energy gain per 
cycle < 2. In the first cycle, the energy gain is large as the particle 
population injected upstream toward the shock is isotropic, hence 
in the first half-cycle u— >d the energy gain ~ r rc i. 

The average energy gain per cycle u— >d— >u is (pf/pi) — 1.7 
to within ±0.1 for the various values of rigidity and shock Lorentz 
factor; this gain tends to diminish with increasing r s h albeit with 
a weak slope. A similar behavior has been observed in the case of 
isotropic downstream turbulence (Lemoine & Pelletier 2003), the 
gain decreasing from ~ 2.0 at F s h = 2 to ~ 1.9 at r s h 3> 1. 

Finally, using the method described in the previous section and 
the probability data collected during the Monte Carlo simulations, 
one can simulate the acceleration process itself and constructs the 
accelerated particle population. The result is presented in Fig. [3] 
This figure reveals that the sub-populations that escape at each cy- 
cle 2n + 1, and whose spectrum is roughly a gaussian centered on 
an energy \/2r s hPo<?™-^d->u (po injection energy) and amplitude 
oc (1 — (Prct)) n , add up to form a featureless power law spectrum 
of index s (at p S> po)- 

The measured spectral index s is shown as a function of shock 
Lorentz factor in Fig. [6] The comparison of these results with those 
obtained for isotropic scattering downstream shows that the inclu- 
sion of shock compression leads to a steeper accelerated spectrum 
at all values of T s h. One can understand this by noting that the 
compressed turbulence leads to lower average return probabilities 
and slightly lower energy gains than those obtained for isotropic 
turbulence (see Lemoine & Pelletier 2003). 

The present results do not settle whether the spectral index 
reaches an asymptote at large shock Lorentz factors, but at the very 
least, as r s h 3> 1 it appears to evolve very weakly close to a value 
s~ 2.6 -2.7. 



3.3 Acceleration timescale 

The present simulations provide a direct measurement of the ac- 
celeration timescale t aC c(e) at energy e, which is defined as the 




Figure 5. Accelerated spectrum of particles escaping downstream times 
momentum squared vs momentum (thick solid line); in thin solid lines, the 
accelerated populations that escape downstream at each cycle. 
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Figure 6. Spectral index vs shock Lorentz factor. 

u^d— >u cycle timescale in the upstream rest frame divided by the 
mean energy gain: 

tu|u(e) + r sh td|d(e/r a h) 



^acc (c) 



(14) 



<?u^d^u 

where t u | u and t&u are the upstream and downstream return 
timescales measured in their respective rest frames. 

The upstream return timescale i„| u ~ 10tL,u/r a h (th,u up- 
stream Larmor time) up to a weak residual dependency on the rigid- 
ity, as shown in Fig.Q A fit that is accurate to a few percent over the 
range of rigidities and for r sh > 5, is: t u | u ~ 14t L ,uP°' 19 / r sh 85 - 
These results agree with and confirm the expectations of Gallant 
& Achterberg (1999) and Achterberg et al. (2001) who argued that 
tulu oc l/r s h since the particles are promptly overtaken by the 
shock when they have been deflected by an angle of order 1 /r s h. 

The downstream return timescale is plotted vs rigidity for var- 
ious shock Lorentz factors in Fig. [8] This figure shows that the 
return timescale t d | d ~ 3 — 4_Rl/c at low rigidities p <C 0.1 
and r s h ^> 1. The uncertainty in the numerical pref actor con- 
tains a weak residual dependence on the shock Lorentz factor 
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Figure 7. Upstream return timescale in units of Larmor time averaged over 
angular distribution vj shock Lorentz factor, for various rigidities as indi- 
cated. The dotted line T rc t = 12T 1 L/r s i 1 is shown as a guide the eye. 
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Figure 8. Downstream return timescale in units of Larmor time averaged 
over angular distribution vs rigidity for various shock Lorentz factors as 
indicated. 



*d|d cx r~" u * (F sh > 1). Note that for the moderately relativis- 
tic shock r s h = 2.2, the return timescale also contains a weak 
dependence on rigidity, t d |d oc p -0 ' 13 approximately, which disap- 
pears at larger shock Lorentz factors. Here as well one can interpret 
the behavior of td|d as the result of reflections of particles on the 
compressed turbulence: the first scattering takes place on a Larmor 
timescale and flings the particles back to the shock with probability 
~ 0.3; if the particle does not return to the shock after this first 
scattering, the probability of doing so at subsequent scatterings be- 
comes negligible for two reasons: the shock moves away at high 
velocity ~ c/3, and the enhancement of the tangential components 
of the magnetic field prohibits efficient transport along the shock 
normal. 

At high rigidities, the downstream return timescale increases; 
the increase is all the more pronounced as the shock Lorentz factor 
is high. This is related to the difference between p« and p; as r s h 
increases, at a fixed value p ~ 0.1, the effective pn becomes larger 
than 0.1 — 1 and the particles leave the range where resonant scat- 
tering with turbulence modes along the shock normal is possible. 



The scattering time thus increases, see Section 4.1, and so does the 
return timescale. 

Finally, the downstream return timescale can be written as 
*d|d(e/r sh ) ^ 4R L , d (e/r sh )/c ~ V3R L , u {e)/r 2 sb , since the 
magnetic field strength is amplified by ^/2/3R. Hence the down- 
stream return timescale measured in the upstream rest frame is 
slightly smaller than the upstream return timescale, and the total 
acceleration time: t acc ~ 10i?L|u/r s h- Interestingly, the accelera- 
tion timescale becomes a fraction of a Larmor time at large shock 
Lorentz factor, which may allow acceleration to an energy limited 
by confinement arguments, in particular Rl < L ma x, rather than 
by energy losses. 



4 DISCUSSION 

4.1 Relaxation length and transport in anisotropic 
turbulence 

The present study has hitherto assumed that the downstream turbu- 
lence is successfully described by the direct compression of the up- 
stream turbulence through the shock jump conditions. However one 
must expect this anisotropic compressed turbulence downstream to 
relax on a timescale r ro i (defined in the downstream rest frame), 
hence on a length scale / rc i = /3 s h|dCr rc i downstream. Particles 
will then experience this compressed turbulence during their jour- 
ney downstream provided the average distance traveled from the 
shock front i tr ~ id|d/2 < l T ei, or 



^d|d i; 2/3 s h|d T "r< 



(15) 



Otherwise the particles reach the point where the turbulence 
anisotropy has relaxed and the previous considerations do not hold; 
however, as we now argue, the above inequality is generally satis- 
fied in relativistic shocks for which the magnetic field is dynami- 
cally unimportant. 

As discussed in Section 4.3, the return timescale £ d | d ~ 
3 — ARl/c for r s h 3> 1 and 2-kRl/ 'Z/ ma x <C 1- There is a resid- 
ual powerlaw dependence on both rigidity and shock Lorentz fac- 
tor but whose power indices are < 0.1, see Section 3.3, which we 
can neglect for now. A simple but somewhat naive estimate for the 
(scale dependent) relaxation timescale is r rc i ~ (kvA)~ , where 
k = 2-k/I is the eddy wavenumber (related to the eddy size I), 
and va is the Alfven velocity. Since particles of Larmor radius 
i?L diffuse through resonant interactions with turbulent modes of 
wavenumber k ~ 1/Rl, the effective relaxation timescale to be 
considered is r ro i ~ Rl/va - The inequality Eq. ll5l is thus satisfied 
when va/c < 0.2 (for T s h 3> 1), which agrees with the hypothesis 
made in Section 2 that the magnetic field is dynamically unimpor- 
tant. Interestingly, one may show that the bound on ?;a is more 
stringent for non-relativistic shocks, since the return timescale (for 
isotropic scattering at least) scales as t d |d c< t S catt//3 s h|d in that 
case. 

The estimate for r re i is likely to be conservative since the eddy 
turn over rate, which gives a refined estimate of the relaxation time 
on a scale k, reads: r t _ ~ (fcffc) -1 , where Vk is now the tur- 
bulent velocity on the scale k, which decreases with increasing 
k; for a Kolmogorov spectrum, Vk oc k~ 5 ^ 3 . With this new esti- 
mate T rc i ~ (fcwA) _1 (fcimax/27r) 5 ^ 3 , the previous condition on 
t rct reads: p < 0.02(/3 sh | d /uA) 3 - Since 7?l ~ I/ max /27r marks 
the maximal energy reached in all likelihood, due to loss of con- 
finement for p > 1, inequality Eq. ^|is valid at all rigidities if 
Va i$ 0.1c. Particles that are accelerated at the shock wave thus 
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Figure 9. Typical trajectory for a particle propagating downstream in shock 
compressed turbulence with shock Lorentz factor r sll = 3.7 (left) and 
r s h = 38 (right). The dotted line indicates the trajectory of the shock 
front; the solid line shows the displacement along the shock normal while 
the dashed line gives the displacement in the plane parallel to the shock 
front. 

do not travel far enough downstream to see anything else than the 
turbulence in its compressed state. 

The transport of particles in strongly compressed turbulence is 
peculiar, as illustrated by Fig.|5| It presents examples of particle tra- 
jectories downstream for two different values of the shock Lorentz 
factor but for the same upstream magnetic configuration; in both 
cases the particle never returns to the shock. The comparison of the 
typical displacement along and perpendicular to the shock normal 
indicates that the particles appear confined in a layer of turbulence 
that lies tangential to the shock plane, and for periods of time ex- 
tending well beyond a Larmor timescale. Note that for r s h = 38 
(right panel), the particle appears to gyrate along a magnetic field 
line located in a plane parallel to the shock front. In this particular 
case, the magnetic field configuration is locally transverse. If the 
trajectory is followed for a sufficiently long period of time, it will 
depart from a simple Larmor gyration law. Moreover the effective 
pulsation is not constant along the trajectory shown in Fig. [5] but 
varies slightly in a random way. It is also apparent in these figures 
that the characteristic pitch angle scattering timescale along the 
shock normal is of order of a Larmor time, while that measured in 
the perpendicular direction is much larger. This demonstrates qual- 
itatively how particles return to the shock in a few Larmor times by 
reflecting on the compressed turbulent modes. 

In order to better characterize the transport of particles in com- 
pressed turbulence, one may seek the diffusion coefficients in the 
various directions. For our purposes, it is more relevant to study the 
time correlation function of the particle velocities, Cy(r): 

Gy(r) = MT)vj(0)>, (16) 

where the average is to be taken on a large number of trajectories, 
and Vi (r) is the velocity of the particle in the direction i at time r. 
The integration of dj (r) over r leads to the diffusion coefficient 
Dij (Candia & Roulet 2004). In the present case, the correlation 
function is however more relevant since particles never actually 
diffuse downstream before returning to the relativistically moving 
shock. The correlation functions along the shock normal, Cm = 
Czz and perpendicular to the shock normal, C± = (d + C yy )/2 
are shown for various values of F B h and p in Fig llOl This figure 
demonstrates that the early time behavior of the parallel correlation 
function is, to a high degree of accuracy, independent of both T s h 
and p, provided r s h 3> 1 and p <C 1. The time behavior of C« can 
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Figure 10. Velocity correlation function C(t) vs. time (in units of Larmor 
time t l ) for various values of rigidity p and shock Lorentz factor r s ^ . The 
lower thick curves correspond to the velocity oriented along the shock nor- 
mal; in thick solid line, p = 3 ■ 10 — 6 and r s h = 3.65; in thick dashed 
line, p = 3 ■ 10~ 6 and T sh = 38.2; in thick dashed-dotted line, p = 10~ 3 
and r s h = 3.65. The upper dotted lines show the correlation function for 
the velocity components perpendicular to the shock normal; at r 15, the 
upper two are p = 3 • 10 — 6 and r s h = 3.65, 38.2, the lower curve is for 
p = 10" 3 andr ah = 3.65. 

be grossly approximated by C\\{t) ~ cos(27rr/tL) exp(— r/rii), 
and tii gives the scattering time along the shock normal. This fit 
is not reproduced on Fig. llOl for the sake of clarity, and because it 
diverges from the measured curves for t/£l > 10. However, one 
can see by eye that the estimated rii w 3iL- 

On the contrary the velocities perpendicular to the shock nor- 
mal do not decorrelate on timescales as short as rii ; the fall-off of 
C± (r) toward zero is observed (but not shown on Fig. llQI on much 
longer timescales than 10 — 20£l, and the decorrelation time rj_ is 
found to depend on p. This is expected insofar as the scattering time 
in isotropic Kolmogorov turbulence tx/£l oc p -2 ^ 3 (Casse et al. 
2002). Indeed, in Fig. llQI the lower dashed curve corresponds to 
the high rigidity value, and it is seen to fall off more rapidly than 
the other two (that correspond to a same rigidity) at r/fe ^ 15. 
One may also note the slight decorrelation in perpendicular veloc- 
ities induced at early times by the decorrelation of velocities along 
the shock normal. 

This figure thus nicely explains the transport properties that 
were indirectly observed in the previous discussion, namely that 
the scattering time along the shock normal, which is the relevant 
quantity for shock acceleration, is of order of three Larmor times, 
and independent of rigidity and shock Lorentz factor to a good ap- 
proximation. 

4.2 Comparison to previous results 

To our knowledge, there is no existing study of Fermi accelera- 
tion in ultra-relativistic shock waves which includes the effect of 
compression of the upstream magnetic field. One may nevertheless 
find interesting points of comparison in various limits with studies 
by Ballard & Heavens (1992), Ostrowski (1993), Bednarz & Os- 
trowski (1998), Kirkef al. (2000) andNiemiec & Ostrowski (2004). 

Ballard & Heavens (1992) were the first to attempt modeling 
of Fermi acceleration in non-relativistic to moderately relativistic 
(r s h < 5) shocks with shock compressed turbulence by the means 
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of Monte Carlo methods. They found a pronounced steepening of 
the spectra index with increasing shock speed and derived the ap- 
proximate formula s ~ 0.75F s h + 1.25. Although our study con- 
firms the increase of s with increasing values of r s h, the precise 
value of s differs significantly from those of Ballard & Heavens 
(1992), all the more so at large shock velocities. One may probably 
attribute this discrepancy to the modest dynamical range (64— 100) 
that was available at the time of the simulations of Ballard & Heav- 
ens (1992). If the dynamic range is not large enough, the particle 
rigidity always lie close to the upper range of resonance and this 
results in steeper spectra due to increased escape probability. The 
subsequent study by Ostrowski (1993) obtained much harder spec- 
tra than Ballard & Heavens (1992) where suitable comparison can 
be made. In particular, for large turbulence amplitude SB / B = 3, 
and shock Lorentz factors T s h = (2.3, 5.0), Ostrowski (1993) ob- 
taineds = (2.0 — 2.2, 2.2 — 2.3), where the range of values bracket 
different values of the mean field inclination with respect to the 
shock normal. Our results for s indicate slightly larger values for s, 
but the agreement is generally better than with Ballard & Heavens 
(1992). 

Kirk et al. (2000) have studied relativistic Fermi acceleration 
using semi-analytical eigenfunction methods; their results confirm 
the canonical value s — 2.23 in the case of isotropic scattering. 
They have attempted to address the effect of anisotropic scattering 
downstream using an analytical description of compressed turbu- 
lence and an analytical estimate of the diffusion coefficient. They 
concluded that anisotropy does not affect significantly s, a result 
which is clearly at odds with the present study. The source of the 
discrepancy lies probably in the modeling of downstream diffusion 
by Kirk et al. (2000): as should be obvious from Figs. 1101 particles 
propagating in strongly compressed turbulence downstream do not 
actually diffuse downstream but rather turn back by reflecting on a 
compressed magnetic layer. 

The most detailed study to date is that of Niemiec & Ostrowski 
(2004), who have studied Fermi acceleration in moderately rela- 
tivistic shocks (r s h < 5) by Monte Carlo integration of the particle 
trajectory in a magnetic field with a large dynamic range. Among 
the results obtained, the authors quote a generic non-power law be- 
havior of the accelerated spectrum: the spectra generally appear 
harder close to the cut-off (p ~ 1) than well below the cut-off. 
Niemiec & Ostrowski (2004) conclude that this effect is probably 
related to the finite dynamic range: close to the cut-off, the prop- 
agation regime in compressed turbulence differs from that in the 
inertial range of resonance, and one indeed finds a different return 
probability or mean energy gain, both of which control the value of 
s. Our Fig.|5|reveals a smooth powerlaw behavior at energies well 
beyond the injection point; this is expected on the grounds that the 
conditional return probability histograms used to model Fermi ac- 
celeration is itself rigidity independent. We have demonstrated in 
the previous Sections that, deep in the inertial range p < 0.1, this 
is a good approximation; however, we have also observed that as 
p tends to larger values, the return probability reveals a slight de- 
pendence on p, and this would make s evolve with p, albeit for 
p > 10 -2 - 10 _1 , had we included this dependence in our calcu- 
lations. 

Our results disagree markedly from those of Niemiec & Os- 
trowski (2004) with regards to the rigidity dependence of the re- 
turn probability in the inertial range. These authors claim to ob- 
serve a pronounced non-monotonic rigidity dependence, which is 
definitely absent from our simulations down to the percent level. 
We note that the noise level in the results of Niemiec & Ostrowski 
(2004) is not indicated in the figures, and might account for part of 



this apparent variability. On theoretical grounds, there is neither ex- 
pectation nor justification for a rigidity dependence of (P re t) on p, 
as long as resonant interactions with the turbulence can occur and 
the shock lifetime is infinite. There is no clear explanation or inter- 
pretation of this observed rigidity dependence of the return prob- 
ability in Niemiec & Ostrowski (2004). Furthermore their mea- 
sured value of (P re t) does not agree with ours: for T s h — 5 and 
SB/B = 3, they find (P ro t) — 0.20, significantly lower than ours. 
This results in a steeper spectrum with s ~ 2.9 ± 0.1, to be com- 
pared with our value s ~ 2.6; one should note that Niemiec & 
Ostrowski (2004) warrant caution with respect to their analysis of 
T s h = 5, as it lies close to the limits of their simulation capabilities. 

We also note that the results of Niemiec & Ostrowski (2004) 
differ signicantly from those of Ostrowski (1993), although the 
method used is similar. For /3 B ^ — 0.5, SB/B = 3 and mean field 
inclination * = 45°, Niemiec & Ostrowski (2004) find s = 2.7 
while Ostrowski (1993) obtains s = 2.0 ± 0.1 for the same val- 
ues (but — 50°). This difference persists at larger inclinations, 
s = 2.8 in the former vs s = 2. ± 0.1 in the latter for * = 75°, 
other values unchanged; it also persists at larger shock velocities, in 
particular for (3 sh = 0.9 (F sh = 2.3) and * = 45 - 50°, s = 2.5 
in the former vs s — 2. ± 0. 1 in the latter. Here as well the source of 
the discrepancy remains unknown. Overall, and where comparison 
can be made, our results lie halfway between those of Ostrowski 
(1993) and those of Niemiec & Ostrowski (2004). 

As stressed at the end of Section 2, the present method offers 
a significant gain in signal when compared to direct Monte Carlo 
methods but it cannot take into account the possible correlations be- 
tween upstream and downstream trajectories due to the correlations 
between upstream and downstream magnetic fields at the point of 
shock-crossing. This remark, when taken together with the above 
comparison to previous work, shows the need for more exhaustive 
studies of relativistic Fermi acceleration with shock compressed 
turbulence, including mean magnetic fields of various strength and 
obliquity, and using both the present method and direct Monte 
Carlo methods. 

4.3 Applications to astrophysical shock waves 

Gamma— ray bursts, with their shock Lorentz factor in excess of 
100 are ideal candidates to test theories on particle acceleration in 
ultra-relativistic flows. In the standard fireball model (see Meszaros 
2002, Piran 2004 for reviews), the prompt 7— emission is inter- 
preted as the product of synchrotron emission of electrons accel- 
erated in the internal shocks with Lorentz factors r s h ~ 2 — 5 in 
the comoving wind frame. The afterglow emission is interpreted as 
the synchrotron light of electrons accelerated in the ultra-relativistic 
shock wave with r s h ~ 300, that itself results from the interaction 
of the 7— ray ejecta with its environment. 

The spectral indices of the shock accelerated electrons derived 
in both cases are s ~ 2.3. This has been interpreted as a dramatic 
confirmation of our understanding of relativistic Fermi accelera- 
tion since it agrees with the canonical value s = 2.2 obtained for 
isotropic downstream scattering. However, as should be clear by 
now, this "agreement" rather reflects our poor understanding of the 
acceleration process: the inclusion of shock compressed turbulence, 
which should be seen as a refinement of the theory, leads to steeper 
spectra, see Section 4, with s ~ 2.4 — 2.6 for r s h = 2 — 5, and 
s ~ 2.6 - 2.7forr sh > 1. 

The difference is not as significant in the case of inter- 
nal shocks than for afterglow observations. As a matter of fact, 
isotropic scattering downstream predicts a value s ~ 2.1 — 2.2 
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for T s h = 2 — 5 (Lemoine & Pelletier 2003), hence it could not 
account reasonably well for the dispersion observed in the spec- 
tral slopes of 7— ray prompt emission. However it is possible that 
the inclusion of a mean magnetic field component with varying in- 
clinations and, possibly varying turbulence level, could reproduce 
this dispersion. Moreover it is not yet established whether the 7 
radiation results from synchrotron emission of shock accelerated 
electrons; other radiating processes (e.g. Piran 2004 and references 
therein) or magnetic reconnection events in the flow (e.g. Lyutikov 
& Blandford 2003) are likely possibilities. 

Concerning the discrepancy of the present spectral index with 
that inferred from afterglow observations, one must note that the 
present study is limited to the case of pure Kolmogorov turbulence 
upstream, which idealizes the limit SB/B 1. However, judging 
by the comparison with Niemiec & Ostrowski (2004), one does not 
expect the inclusion of a coherent component to help resolve this 
discrepancy, as these authors have observed a steepening of the ac- 
celerated spectrum with decreasing turbulence level SB/B. In the 
ultra-relativistic regime, this trend should be exacerbated, since as 
SB/B decreases, one approaches the perpendicular shock accel- 
eration limit where Fermi acceleration becomes inefficient (Begel- 
man & Kirk 1990). Our simulations should thus provide a conser- 
vative lower bound to s for the case that includes a mean magnetic 
field. 

The discrepancy might be attributed to the nature of the tur- 
bulence, in particular to the assumption of Kolmogorov turbulence. 
Again, the work of Niemiec & Ostrowski (2004) suggests that the 
turbulence spectral index has an effect on s, although there is not 
enough simulation data to pinpoint what the exact correlation is. 
These remarks indicate the need for more exhaustive studies that 
investigate various turbulence spectra. Interestingly, this suggests 
that prompt and afterglow observations of 7— ray bursts might be 
giving us information on the properties of the turbulence behind 
the shock front; the upstream turbulence does not play any role in 
the ultra-relativistic limit as a result of the limited amount of time 
that a particle spends upstream before being overtaken by the shock 
front. 

The interpretation of 7— ray bursts afterglows as synchrotron 
emission by shock accelerated electrons requires that the magnetic 
field intensity at the shock front be significantly higher than the 
average interstellar value (e.g. Piran 2004 and references therein). 
The nature of downstream turbulence would then be directly related 
to the amplification process, which might manifests itself indirectly 
in the spectral slope. The proposal of magnetic field amplification 
by the two-stream Weibel instability (Medvedev & Loeb 1999) has 
recently triggered a lot of interest. This instability seems able of 
explaining the high value of B required, although debate on the 
subject is not closed, see Wiersma & Achterberg (2004). In any 
case, the Weibel instability amplifies the magnetic field in the trans- 
verse plane to the shock normal and on very small spatial scales 
(Medvedev & Loeb 1999) . If the magnetic field on scales larger 
than or equal to the Larmor radius of a typical accelerated particle is 
thus not amplified by the instability, one should expect acceleration 
to proceed as presented here, and the discrepancy should remain. 
These considerations may suggest that the amplified magnetic field 
structure differs from that proposed by Medvedev & Loeb (1999) 
or that the amplification mechanism itself is different. In this re- 
gard, we note that recent numerical studies on the Weibel instabil- 
ity suggests that powerlaw acceleration may occur independently 
of the Fermi mechanism due to the presence of electromagnetic 
currents downstream (Hededal et al. 2004). It should also be noted 
that stochastic acceleration in the downstream turbulence, which 



has not been accounted for in the present study, could play a sig- 
nificant role in reshaping the accelerated spectra, as suggested by 
Virtanen & Vainio (2005). 

On a different line of thought, one should point out that the 
value of the spectral index derived here turns out to be in very 
good agreement with that required to fit the ultra-high energy part 
of the cosmic ray spectrum at energies E > 10 18 eV, namely 
s ~ 2.6 — 2.7, when one assumes that the sources are distributed at 
cosmological distances and do not evolve too strongly with redshift 
relatively to the cosmic star formation rate (Berezinsky et al. 2002, 
Berezinsky et al. 2005, Lemoine 2005). 



5 SUMMARY 

We have conducted a study of Fermi acceleration at relativistic and 
ultra-relativistic shock waves, considering the effect of the shock 
compression on the downstream magnetic turbulence. The numeri- 
cal simulations are based on Monte Carlo methods of particle prop- 
agation in realistic magnetic fields described by sums of plane wave 
modes. The numerical technique differs from the standard Monte 
Carlo modeling of Fermi acceleration in that it measures the rele- 
vant statistical laws of particle transport on either side of the shock, 
and uses these probability laws together with the Lorentz transform 
from one frame to the other to reconstruct the acceleration process. 

The turbulence was assumed to be described by pure Kol- 
mogorov turbulence upstream, a situation which idealizes the limit 
SB /B ^> 1. The main effect of the compression with respect to the 
case of isotropic scattering is to steepen the accelerated spectrum to 
a slope s ~ 2.6— 2.7 in the limit r s h 3> 1, as a result of a decreased 
return probability. This latter effect is induced by the compression, 
which amplifies the magnetic field in the transverse direction to the 
shock normal: particles that enter downstream are trapped on a Lar- 
mor timescale in a compressed turbulence layer and cannot recross 
the shock unless they turn back within a few Larmor times. Con- 
sequently, the acceleration timescale is dominated by the usptream 
residence time, and can be as short as i acc ~ 10£l /T s h (upstream 
frame). We have also argued that the accelerated particles do not 
travel far enough downstream before returning to the shock to ex- 
perience a turbulence that has relaxed to near isotropicity. 

The derived slope does not agree with that inferred from ob- 
servations of 7— ray bursts afterglows, which indicate s ~ 2.3. This 
inferred value is generally accepted as a success of Fermi acceler- 
ation in the relativistic regime whose predicted canonical value is 
s ~ 2.2 — 2.3. However, this result only holds for isotropic scat- 
tering downstream, whereas the inclusion of realistic shock jump 
conditions, as done here, makes downstream turbulence strongly 
anisotropic and the spectra markedly steeper. The resolution of this 
discrepancy may be tied to the necessary but unknown amplifica- 
tion mechanism of the upstream magnetic field. 
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